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Abstract 

Quantum state preparation through external control is fundamental to established methods in 
quantum information processing and in studies of dynamics. In this respect, excitons in semicon- 
ductor quantum dots (QDs) are of particular interest since their coupling to light allows them to 
be driven into a specified state using the coherent interaction with a tuned optical field such as 
an external laser pulse. We propose a protocol, based on adiabatic rapid passage, for the creation 
of entangled states in an ensemble of pairwise coupled two-level systems, such as an ensemble of 
QD molecules. We show by quantitative analysis using realistic parameters for semiconductor QDs 
that this method is feasible where other approaches are unavailable. Furthermore, this scheme 
can be generically transferred to some other physical systems including circuit QED, nuclear and 
electron spins in solid-state environments, and photonic coupled cavities. 
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Whether one is interested in methods of quantum optics or computation, or more general 
many-body dynamics, it is important to be able to prepare interacting few-level systems in 
initial states from pure product states to more complex entanglements. Excitons in low- 
dimensional semiconductors exemplify such a level structure. They are of interest both 
because their strong coupling to light provides communication with the external world, 
and because photon-mediated interactions between excitons can be substantial. However, 
the lack of degeneracy due to the (accidental or deliberate) inhomogeneous spread of the 
energy levels becomes a challenge for addressing different levels - either independently, 
for example to invert a single one, or in pairs, for entangling a specified set. Control of 
the population of excitonic states in quantum dots (QDs) has been demonstrated using 
resonant transform-limited light fields such as vr Rabi pulses P-H]. Although effective in 
certain cases, Rabi flopping has the disadvantage for population transfer that the final state 
is sensitive to variations in dipole couplings, dot energy and disorder. Such disorder becomes 
a serious challenge if one wishes to scale this approach beyond a single quantum dot, and 
similar limitations arise in many other systems. More recently, population transfer has been 
achieved [HI E] using a protocol called adiabatic rapid passage (ARP) [7] . Within this scheme 
the frequency of the laser pulse is swept through that of the target state and, if whole process 
is performed adiabatically, this target state will be populated with high efficiency. Unlike 
the Rabi approach, ARP is largely unaffected by all the variations previously mentioned, 
and this has prompted theoretical proposals for its use in preparing populations in ensembles 
of QDs [S], as well as proposals to realize a non-equilibrium Bose-Einstein condensate p-[TU] 
and quantum operations [TlHTT] . 

The aim of this paper is to show how ARP may be generalized to create entangled states 
in a disordered system, comprising a set of two-level systems coupled together in pairs. Such 
an ensemble could be realized using quantum dot molecules, with a possible intra-molecular 
coupling mechanism being resonant Forster transfer [131 HHHSD]- Implementations would 
also be feasible in other systems such as photonic coupled cavities with photons tunneling 
between neighboring cavities and having onsite interaction [21] , in arrays of superconducting 
qubits with an exchange coupling mediated by virtual photons [221 123], or using electron 
spin states in coupled semiconductor impurities [21|. 

The important feature of our proposal is that it requires neither precise control, nor 
reproducibility, of the pairs. Instead we exploit the disorder to allow those pairs with 
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entangled states to be spectrally identified within an ensemble. The robustness of ARP 
then allows a pulse to be constructed which drives those pairs into their entangled states, 
without creating other excitations. Thus, as shown in Fig. [2| the entanglement of formation 
per excitation is very close to one, even in a strongly disordered system. Our approach 
represents a significant simplification and improvement of the protocol required for the 
production of entangled states in realistic systems, since previous proposals have relied on 
coupling two states through a further level |15], or consider only a single and fine-tuned 
system [nHISl ESI • 

To simplify the notation we consider the limit in which the states of the individual 
quantum dots form two-level systems that may be represented in terms of a Bloch vector or 
pseudospin: each dot either contains no exciton (spin down) or a single exciton (spin up). 
The generalization to allow for multiple states on the dot, for example due to the exciton 
spin structure, is straightforward. The relevant Hamiltonian for two coupled dots, in a frame 
rotating with the time- dependent frequency uj{t) of the ARP pulse and in the rotating wave 
approximation, is then {h = 1): 



i=l,2 
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where ii = Ei — uj{t), Si is the transition energy for dot i = 1,2, g(t) the amplitude of the 
driving field used to perform ARP, cij the Pauli operators describing the state of dot i, and jx 
an exchange interaction between the two sites {h = 1). The ARP pulse can be decomposed 
into its amplitude and frequency, i.e., F{t) = g{t)exp (^i f uj{t')dt'). For definiteness we 
consider a Gaussian linearly chirped pulse: 

F{t) = (/(t)e'"(*)* = (?oe-*'/'"'e*"(*)* , oo{t) = ooo + at , (2) 

with go the pulse amplitude, r the temporal width of the pulse, uq its central frequency, 
and a the chirp. We take the pulse width to define units of time and energy, and use 
dimensionless parameters ar^, g^r and Jtt. 

In general, ARP schemes are generalizations of the Landau-Zener problem j271 |2B] to a 
time-dependent and controllable mixing g{t) between the levels. Consider first noninteract- 
ing dots (jT = 0): if the driving frequency sweeps through the transition frequency of a dot, 
Ei passes through zero, and the lowest energy eigenstate changes from an empty site to being 
fully occupied. In the Landau-Zener problem there is a constant hybridization go between 



these levels to create an anticrossing, and a system in the initial ground state will be driven 
into the excited state with probability Pinv = 1 — exp(— (/q/q;). Thus as the adiabaticity is 
increased, by reducing the chirp a or increasing the hybridization g^^ the final population 
rapidly approaches one. In the case of ARP, the coupling term gii) is not constant and, in 
order for the transition to occur, the two levels of the system need to be coupled long enough 
so that the character of the eigenstates changes smoothly, i.e., adiabatically. This gives the 
criterion ar^ ^ 1 Furthermore, the ARP scheme can be made spectrally selective. One 
can choose to excite either one dot or the other by timing the center of the optical pulse to 
occur when the chirp has driven a particular level to zero, and ensuring that the pulse is 
negligible at the other crossing. The available discrimination in the energy levels is thus of 
order ^/r. 

This discrimination could be used to create entangled states in the ideal case of a pair 
of identical two-level systems [TTHT3]. In this case the Hamiltonian is, in the ab- 
sence of the driving field, diagonal in the singlet /triplet basis {15*), |T_), |To), with 
\S) = il/V2) (I ;t) - I U)), \T^) = I U), \To) = il/V2) (I ;t) + I U)), \T+) = I tt). The 
intermediate states \Tq) and 15) between the ground state and the fully occupied state |T+) 
are spatially entangled if the coupling dominates over the detuning, Ae = £i — £2 ~ jr- For 
a single pair of identical dots (^i = 62 = Sq), the use of ARP to create this entangled state 
is illustrated in Fig. 1(a) which shows the time evolution of the eigenvalues of the Hamil- 
tonian Eq. (1) when the levels are coupled by a linearly chirped Gaussian pump {qot = 5, 
colored continuous lines) and uncoupled [qot = 0, thin dashed lines). By chirping through 
the level crossing (A) between the states |T_) and |To) (at uq = Eq — jr) with a pulse dura- 
tion such that coupling switches off before the \Tq) — |T+) crossing (B), the system will be 
adiabatically driven and left in the entangled triplet state \Tq). 

We now consider the parameters required for this procedure for a single degenerate pair. 
In this case the population of the two-dot system in the desired final state is equally shared 
between | and | ^t); so p^i = = 0.5. The element of the density matrix p^^ (= p^^), 
can thus be used as a measure of the insensitivity of the entanglement to variation of inter-dot 
coupling and chirp. This is illustrated in Fig. 1(b) which shows p^i, obtained by numerical 
solution of the Liouville equation dp/dt = —i[H,p], as a function of the dimensionless 
parameters ar^ and Jtt. The region where p^i = 0.5 is seen to extend over a large range of 
values of chirp and exchange coupling. The required parameters are fully compatible with a 
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FIG. 1. (Color online) (a) Time-dependent energy levels of a coupled system of two identical dots. 
The continuous colored lines indicate the energy levels with a Gaussian driving pulse {g^r = 5) of 
duration r, while the thin dashed lines are the undriven levels (170 = 0). The chirp and interdot 
couplings are ar^ = 2 and jrpT = 4. The central frequency of the ARP pulse is Eq — Jt, resonant 
with the transition between the ground state |r_) and the entangled state {Tq) (see point A), (b) 
Final density matrix element p-^^, calculated as a function of the dimensionless linear chirp ar^ 
and exchange coupling Jtt. The white region (p-^^ = 0.5) shows the range of values of chirp and 
exchange coupling where the pair is driven into the entangled state \Tq) with high probability. 
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linearly chirped pulse similar to that used to invert a single semiconductor QD [5]: in that 
case we have considered a transform-limited pulse width of 2 ps and a chirped temporal 
width r = 4.5 ps. Thus for ar^ ^ 2 we require a dimensionless coupling j^r > 4 [see 
Fig. 1] for the ARP transition to the entangled state to occur, which corresponds to a 
value of Forster coupling of 0.6 meV. Previous studies have estimated an upper limit of 
10 meV for Forster coupling in semiconductor quantum dots [20]. Thus the scheme could 
be implemented using two stacked (vertical) QDs at a distance of few nanometers, coupled 
by Forster energy transfer, but with no single-particle tunneling. Such conditions can be 
achieved in InAs/GaAs coupled QDs [221 ED]- The scheme would also apply with single- 
particle tunneling or indirect excitons, provided only that target entangled states can be 
identified. 

We now show how this adiabatic protocol can be generalized to generate entanglement 
in realistic non-degenerate systems such as ensembles of QDs. In this case the entanglement 
pulse can be spectrally tuned to address a specific pair of dots, and several pairs can be 
entangled within an ensemble by superposing such pulses. 

It is important to recognize that the requirement of exact degeneracy of the uncoupled 
transition is relaxed up to the magnitude of the coupling energy. This affords a route 
to practical realizations of the scheme, as the coupling energy and level splitting can be 
traded to optimize the probability of producing an appropriate double-dot structure. In the 
following we consider an ensemble of such systems, modeled here by an average coupling 
strength and having an inhomogeneous distribution of energies. In each coupled pair the 
energy levels are not degenerate, and only pairs with detunings smaller than the interdot 
coupling strength can be entangled; these specific pairs must be identified spectroscopically 
prior to ARP manipulation. Spectrally-selected components of a broad-band pulse, each 
close to resonance with a particular chosen pair, can then be chirped in the same linear 
optical process, such as a grating-based delay stage. We test this scheme simulating an 
ensemble consisting of thirty pairs of dots with energies taken from a Gaussian distribution 
of standard deviation 10 meV, and coupled by an average coupling strength Jt =2 meV. In 
one typical realization taken as an example, three couples can be entangled as their energy 
levels are detuned by an amount smaller than jT/2. To evaluate the entanglement produced, 
we use the entanglement of formation (EOF), which is a widely accepted measure of the 
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entanglement for bipartite states. For pure states it is the entropy [31] : 

Em) = S{p), (3) 

where S{p) = — Tr (^p^'^'^log2p^°'^) is the von Neumann entropy and p^'^'^ = Tr^^^\"^) {'^\ is the 
reduced density matrix obtained by tracing the whole system density matrix |\1')(^| over one 
of the two subsystems of which the pure state |\E') consists. In this case the two subsystems 
are the two paired dots. The EOF (or entropy) ranges from zero (for a product state) to 
log2A^ for a maximally entangled state of two A^-state particles. Hence, the EOF of the 
triplet state |To) = (l/-\/2)(| ti) + I it)) (which is a maximally entangled Bell state) is equal 
to 1. Furthermore, the properties of the EOF as a natural measure of the entanglement 
include that the entanglement of independent systems, such as pairs of coupled dots in an 
ensemble, is additive. 

Figure 2 shows the EOF and the total excitation in our example realization calculated as 
a function of the strength of the applied chirped Gaussian pulse. The excitation induced in 
the system after the application of an external pulse has been evaluated as X = Tr(pX), the 
excitation operator X given by X = (cr^ + l)/2 + (cr2 + l)/2 and (j| 2 being the two-particles 
spin operators for dot 1 and 2. We have used r = 15 ps and a = 0.01 ps^^ (a;r^=2.25) 
yielding a small value of ar - which determines the energy range spanned by each component 
of the ARP pulse - while keeping the whole process in the adiabatic regime ar^ > 1, so that 
few states apart from the entangled ones will be excited. The total entanglement (dashed red 
curve) always deviates from the ideal value of 3, i.e., the total number of triplet states, due to 
the unavoidable excitation of other states which are not entangled, but yield a small amount 
of entropy. Furthermore, as the intensity of the pulse increases, other level crossings become 
adiabatic and other states such as the full inversion | tt) can be populated, thus explaining 
the behavior of the excitation (continuous black) curve. However, the EOF normalized with 
respect to the total excitation (dashed-dotted blue curve), which is an effective measure of 
the total entanglement produced in the system, is close to the ideal value of 1 (~ 0.85 for all 
the different pump strengths considered). In order to validate these results we have simulated 
50 different realizations of Gaussian-distributed coupled dots, and found an average value 
for the normalized EOF always ~ 0.8 for several realistic values of pulse strengths. 

In conclusion, we have shown the benefits of ARP for the generation of entangled stated 
in ensembles of coupled quantum systems such as inhomogeneously distributed QD excitons. 
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FIG. 2. (Color online) Predicted entanglement generated when an ensemble of 30 coupled pairs of 
dots is driven by an ARP pulse. The pulse is constructed to drive the 3 strongly-coupled pairs in the 
ensemble into their entangled state (see text). The red dashed curve shows the total entanglement 
entropy (or entropy of formation), the black solid curve the total excitation, and the blue dot- 
dashed curve the entanglement entropy normalized by the total excitation. The horizontal axis is 
the strength of the driving field, defined as twice the peak Rabi splitting in a single dot, given in 
units of 1/r {h = 1). 



Our calculations, based on realistic values for Forster coupling and dot distribution, provide 
a feasible route for the realization of entanglement in solid state systems, which is more 
practical than other approaches due to the flexible and tunable parameters of ARP. In 
principle, the same methodology can apply to microcavity Josephson junctions, to coupled 
electron spins in semiconductors (using ESR), and to coupled nuclear spins (via NMR). 
The latter arena is the genesis of the ARP technique, though we are not aware of it being 
explicitly used in the form proposed here. 

C.C. and R.T.P. acknowledge support from EPSRC under Grant EP/F040075/1. P.R.E. 
acknowledges support from Science Foundation Ireland under Grant 09/SIRG/I1592. P.B.L. 
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